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Abstract 

We investigate the effects of given pressure gradients on hydrodynamic flow equations. We obtain results 
in terms of implicit solutions and also in the framework of an extra-dimensional formalism involving the diffu- 
sion/Schrodinger equation. 
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£5 ' 1 Introduction 

O ' 

CO ' There has recently been renewed interest in the similarities between the equations of fluid dynamics and those of 
quantum mechanics PP-0- This is a venerable subject UJ. In a previous article |1U| we have shown how any solution 
of the multi-variable, multi-dimensional generalization of the simplest flow equation, the Euler-Monge equation, is 
related to a particular type of solution to a linear diffusion or heat equation, in twice the number of "spatial" 
dimensions. But, as is well-known, the heat equation when complexified is just the Schrodinger equation. So 
the method of may be viewed as a new transformation relating the Euler-Monge and dimensionally-doubled 
C* { ■ Schrodinger equations, which differs from the transformation of Madelung. 

Here, we extend this result, especially in the case of one space dimension, to a flow equation driven by some 
particular, given pressure terms (i.e. so-called "body forces"). When the given pressure-gradient is either t or 
x independent, we obtain the general solution for such pressure driven flows in closed but implicit form, and we 
relate these solutions to those of driven diffusion/ Schrodinger equations. The corresponding Schrodinger equation 
contains a potential whose form is related to the pressure-gradient. We also find explicit series solutions in a variety 
of situations. Special attention is given to constant and linear pressure-gradients. We discuss the complexification 
required to obtain the Schrodinger equation, and we emphasize analogies between the heat equation and the Wigner- 
Moyal equation for the evolution of quantum mechanical densities. In addition, we show how solutions of a related 
Bateman equation may be constructed. 

Should there be any lack of motivation to consider such matters, perhaps a brief historical overview is warranted. 
The Euler-Monge (E-M) equations embodied in fluid dynamics first appeared in 18th and 19th century studies of 
both fluid flow |llj and analytic geometry |12| . Riemann took up a study of the equations in the context of gas 
dynamics, discussing the equations as a theory of invariants \YA\ (for a modern textbook treatment, see |14j). His 
approach is widely applicable to almost all nonlinear flow problems, although it does not triumph over turbulence. 
A systematic modern discussion of the E-M equations that synthesizes ideas from both geometry and invariance 
theory can be found in the review by Dubrovin and Novikov |15|. Most contemporary texts and reviews stress the 
universal role played by these nonlinear transport equations in accordance with Whitham's theory |16| : Essentially 
all nonlinear waves, even those in dispersive and dissipative media, involve E-M equations, or simple variants of 
them, if the nonlinear wavetrains are slowly varying. This makes the equations particularly useful for analyzing the 
asymptotic behavior of nonlinear solutions. The E-M equations and their conservation laws also serve as a useful 
starting point in Polyakov's study of turbulence [Q but without yet leading to a general solution of the Navier-Stokes 
equations. In short, the Euler-Monge equations appear widespread across a very broad landscape of physics and 
applied mathematics problems, and therefore it is important to understand their solutions at as many levels as 
possible. To that end we shall map some solutions of the E-M equations with given pressure-gradients onto linear 
diffusion/ Schrodinger equations. 
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Suppose we ask for the condition that V(x,t,a) = ^expau(x,t) is a solution of the following linear diffusion 
equation in two (x and a) dimensions |17| with a given "potential" term a x g(x,t) : 

This holds if and only if a nonlinear equation is satisfied by u = - In (aU) : 

du(x,t) , . du(x.t) , . 

We will call this the "pressure driven Euler-Monge equation." The function g is the pressure-gradient, g(x,t) = 
d x p(x,t), as it would appear in the Navier-Stokes equation for an incompressible fluid without viscosity. Here, 
we shall concentrate on pressure-gradients which are given functions, especially those which are time- independent, 
g{x) = d x p (x), and not have p determined by an equation of state. (Hence, our discussion is not truly a description 
of ideal non- viscous fluid flow.) 

The linear equation can always be attacked, at least formally, by the propagator method. However, in 
general, it is not possible to find explicit, closed-form solutions, either to that linear equation or to the associated 
nonlinear equation Pjl. Nevertheless, general implicit solutions to the nonlinear equation can be obtained in closed 
form, as shown below, when the pressure-gradient is either a function of x alone, or of t alone. Moreover, in cases 
where the pressure takes simple forms, special explicit solutions exist, in terms of which the general solution may be 
constructed. 



2 Simple examples 

In this section we discuss some simple cases before treating a general pressure term. The technique we use is a 
variant of the well-known method of characteristics. 

First, consider the case of the constant pressure equation, with g (x) = 0. 

du du 

m= u dx-- (3) 

Two elementary special solutions are given by 

x + ut = , u = . (4) 

Now, to obtain the general solution, just set one of these expressions to be an arbitrary differentiable function of the 
other. 

x + ut = G{u) , (5) 

for any differentiable G, where in (J5J, u = u(x,t). This result is well-known By taking x and t partial 

derivatives, the reader will easily verify that as a consequence of |J5J, u(x,t) must satisfy the Euler-Monge equation 
©■ 

Next, consider the equation resulting from a linear pressure, p (x) = p (0) + kx, hence constant pressure- gradient, 
9 (x) = k. 

du du 



Two special solutions are given by 
Then the general solution is given by 



m= u d~x + k - (6) 

x + ut- ^kt 2 = , u-kt = 0. (7) 
x + ut-^kt 2 = G(u-kt) , (8) 



where again G is any differentiable function, and where u = u (x,t) in this last equation. 
In the case of a quadratic pressure the equation becomes 



du du , 2 , n . 

m= u dx + kx > (9) 
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with a linear forcing term k 2 x. A particular solution is now given by 

kx cos(fct) + u sm(kt) = . (10) 

(Dividing by k and taking the limit k — > recovers the first of Q.) Now we can construct another solution by 
translating kt by -| to give 

fcxsin(fci) — ucos(kt) — 0. (11) 
(In the k — > limit this gives the trivial solution in Q.) The general solution is now given implicitly by 

kx cos(kt) + u(x,t) sin(fei) = G(kx sin(fct) — u (x, t) cos(fct)) , (12) 

with G an arbitrary differentiable function. 

3 Recipe for a general implicit solution 

One way of proceeding with equation J5J, in the case of arbitrary time- independent g (x) , is not to consider u as a 
function of x and t, but to think of t as a function of x and u, i.e. the hodograph method [16j . Then we have 

du dt(x, u)/dx du 1 

dx dt(x 1 u)/du' dt dt(x,u)/du 

This transforms equation J3J into a linear equation for t. 

d d 

— (ut(x,u)+x)=g(x)—t(x,u). (14) 

The solution for general g{x) is easily constructed from this. We turn this equation into an integro-differential 
equation by integration with respect to x. 

x 1 d [ x 

t(x,u) = F(v.) 1" — 5- / t(z,u) g(z)dz , (15) 



u u du J 

where F(u) is an arbitrary function of integration. This equation is then formally solved by iteration: 

/ x n / \ x 1 a r x , . . /„. . z x \ 1 a r , , . 1 a r 21 , . , , 

t{x,u)=F(u) h — / dzig(zi) \F(u) + -— / dzig(zi)-— dz 2 g(z 2 )t (z 2 ,u) 

u u ou J \ u / u ou J u ou J 

= F{u)-- + ■■■+ / dzig(zi)-.- dz n9 (z n )(-f) (*»_£»)+.... (16) 

zt J Jo \uduj V tt / 

Since every term in the expansion is a product of a function of x with a function of u, all these terms are easily 
disentangled. The terms linear in F can be evaluated separately to give 

id r , , , r* 1 , , , / 1 9 N 2 



1+ / dz l g{zi)-— + / dz 1 g(zi) / dz 2 g{z 2 ) [-g^l H ) F ( M ) 



- fa (jf **«>)" (il:)" F <"> ^ [,; "" liif ("i - f (-^+IT*»w) 

The (n + l)th term in the remaining series of terms is given by (— l) n+1 (2n — l)!!/it 2n+1 times 



(17) 



Zl 



JO JO 



(n-l)Uo 



dz 1 g{z{) I dz 2 g(z 2 )--- I dz n g(z n )z n = - — ^— I (p(x)-p(z)) n 1 g (z) z dz , (18) 



where the pressure is p(x) = p(0) + dzg(z). These terms can be summed. After an integration by parts, the 
final result is 

t (x, u) = F + A(p(x)-p(0))) - ± f ^ . (19) 

This is the general solution to (fT^|) with boundary value t (x = 0, u) = F (u). 
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The result (|19fl is easily checked to be a solution to (|14f) . We recognize this as again being of the form 
where an arbitrary function of one particular solution of J5J) is set equal to a second particular solution. In 

this case F (^u^Jl + ^ (p (x) — p (0))^ = G{^u 2 +p{x)) involves an arbitrary function of the combination that 

gives a particular time-independent solution, ^u 2 +p(x) — constant, of the time-independent special case of (0, 
hd x (it 2 ) = ~d x p(x), and the remaining terms in t(x,u) implicitly furnish another particular solution of the fully 
time-dependent 

Similar statements apply to the situation where the driving term is a function of only t, and not x. Thus 
du/dt — udu/dx + kit) is solved implicitly by x + ut — J * k (r) tcLt = G (u — k (r) drj , where the latter function 

has as its argument the same form as the particular x-independent solution u — f k (r) dr — constant. 

4 Series solutions 

Building up to the extra dimension Of course, the existence of an implicit solution, while providing a pleasing 
theoretical establishment of integrability, is not much use in many practical applications. So we also look for power 
series expansions in t. Consider hrst the case of a constant pressure-gradient, 

du du , 

Tt= u Tx + k - (20) 

Further time differentiation removes the explicit fc, just as it would for any time-independent g (x), 

d 2 u ldd 2 d.du d,du. 

W = 2did^ U = di {u d-x ] '= Tx^ " (21) 

If the power series is denoted by 

oo 

u(x,t) = J2 tn M%) , ( 22 ) 

n=0 

we have u (x, 0) = u (x) = Uo, u% = k + \-§^u (x) 2 , and so on. That is, 

Id n 

(n + 1)m„+i = -— 22 u J u n-j ' for n > 1 . (23) 

X 3=0 

The resulting expression for d n u/dt n \ t=Q = n\u n for n > 2 is therefore given, for constant k, by 

L«/2J 



§ \2j jl(n-2j)\dx n -i- 1 \ U{x)n 3 ~dx~) • (24) 



Here, the floor function, [ n /2j, is the greatest integer less than or equal to n/2. Hence, 

oo L"/2J 



itx, t) = u(x) + kt+y^ t n —7 — - — -T7 - ^ — -u (x) 

y 1 ' w 4r! j!(n-2j + l)! V2/ <9x"-J v ' 



n-2j + l 



n=l j=0 

L !!/2 J 1 / I. * \ J / »2 \ «-2j 



i { k d Y { d 2 \ 

= + * + i!(n-2i + l)! UtoJ 



it (x) 



(25) 



a=0 



This last expression involving the exponential, evaluated at a = after all differentiations, evokes the extra- 
dimensional approach. 

To follow through on that idea, recall that Hermite polynomials are given by (^Hjj 22.3.10) 

L"/ 2 J 1 L«/2J , , , v j 
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So, taking d/dx of both sides of l|24[l . we find, formally, 

L«/2J 



at™ dx "lt=o / , 

3=0 



(fc/2)' f..n-2j 
{j)\(n-2j}\ dx"-J V U Ux 



[n/2\ 

(i)!(n-2j)! V ~ 77 



fe _9_y / a 2 \ 

2dx) \dxdaj 



n-2j 



a=0 



n L«/2J 



(j)!(n-2j)! 



i-iy 



k d 



Ik 9 
2 THc 



n-2j 



a=0 



\ 2 dx J n\ " n 

This result allows us to perform the sum, to obtain 



/ k 9 
2 Tlx 



(27) 



a=0 



71 = 



exp 



n=0 



2t /fc a 



l -' 



dxOa 



-1 V 2 9i o / fc a 



t /fc JL 



/ fc a 
2 a 



It (x) 



exp 



a=0 



a=0 



a=0 



., 9 e a«(x) _ 1 



dx 



(28) 



through use of the generating function ([THJ, 22.9.17) e 2 ^" 8 = XT=o ^"^n ( z ) ■ Thus 

we conclude 

u(x,t) = kt+e t ^ + ^ kt ^ =kt+e t^- a l( e au( x+ ikt 2 ) _{\ (29) 

V a ) a=0 a V / q=0 

The RHS clearly involves just the time evolution of U, evaluated at the extra dimension boundary, a — 0, but under 
the action of a slightly modified extension of the time-dependent kernel in |10|. such that the initial data has a 
supplementary shift x — > x + ^kt 2 . 

Comparing this to the undriven situation, we see the physical interpretation is simply that of a constant external 
acceleration imposed uniformly on the nonlinear self-interaction of the fluid. If we "fall" along with the self- 
interacting fluid, constantly accelerating with it, the effects of the external acceleration are removed. That is to say, 
from Galilean covariance, if u (x, t) is any solution of Q, then u (x + ^kt 2 , i) + kt is a solution of ©. 

Working down from the extra dimension The previous implicit solutions suggest another way to obtain results 
in terms of an extra-dimensional equation: Incorporate one of the particular solutions into the exponential along 
with u (x, t). For example, define 



(30) 



,( X ,t,a) = i ^Mx,t)-akxt a n(kt) _ {j 

Then we have 

(I " 5§£ ^ k tan kt + + a £)) U (*> *. a ) = e au ^~ ak * tan < fct ) (d t u (x, t)-u (x, t) d x u (x, t) - k 2 x) . (31) 



This gives an equivalence between the Euler-Monge equation with linear pressure-gradients, and solutions, of 
the form given in (|30|l . to the following modified extra-dimensional linear equation. 



d d 2 



d d 

il (x, t, a) = (k tan kt) I 1 -I- x— — ha— ) il(x,t,a) 

ox oa 



s dt dadx ) 

The formal solution to this linear equation is given by the propagator method. 

il (x, t, a) = K (i) il (x, t = 0, a) = K (t) ± (e au{x) - l) 



where the operator kernel obeys 



d 



(32) 



(33) 



(34) 
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That is, the kernel is a time-ordered exponential [T5| given by the usual series of nested integrals 

K W = 1 + E j ( dTl J ldT2 '"J ' 1 dT > (af^ + A;tanfcri ( 1 + x ]L + a ]L 

x (^ + fctanfcT2 ( 1 + ^ + a ^))---(9^ + A;tanfc ^( 1 + x ^ + a ^)) ' ( 35 ) 

This timc-ordcrcd exponential can be worked out as products of ordinary exponentials, through use of the underlying 
algebra. 

[A, C] = = [B, C] , [A, B] = A , (36) 
where A — g ° ga , B — a-^ , C = 1 + — a^. Note that C is central. Writing 

<9 2 f d d\ d 2 ( d\ (3d 
+ (k tan kt) 1 + a;— + a— = + 2 (fc tan fci) a— + (k tan fei) I 1 + a;- a- 



dxda \ dx da J dxda \ da J \ dx da , 

= A + 2 (k tan kt) B + (k tan kt) C , (37) 

we see that we need to solve an evolution equation of the form 

±K(t) = (a(t)A + {3(t)B + 7 (t)C)K(t) , (38) 

with K (t = 0) = 1, where A,B,C obey (|36|l . and where the coefficients a, /3, 7 commute with everything. The 
solution is given by the Baker-Campbell-Hausdorff technique. 

K(i) =exp^B J f3 (r) drj exp a (t) efo > 3 ( T ') dT ' dr^j exp (^C J (t) dr^j 

= exp J^ a ( T ) e ~ J * P ^ dT 'dT^j exp jf /3 (r) dr^j exp ( C J q 1 ( t ) dr ) ■ ( 39 ) 

For the case at hand, a (i) = 1, /3 (t) = 2kta,nkt = 2j (t), f Q /3 (t) dr = — In cos 2 kt — 2 f Q 7 (r) dr, and 
$la{r)e-tt fi ( T ') dT 'dT= ^sin2fct. So 

K (t) = exp ^^-^Isin (2kt)j exp (-Bin cos 2 kt) exp ^-^C In cos 2 fci 

= c^ eX K^ Sin(2fc ^) eXP (^ ln(COS2H) (^ + ^)) • (40) 

The right-most exponential operator rescales x and a in the initial data, and the left-most operator then evolves 
these data as in the undriven case, except for the replacement t — » 4j- sin (2kt). 
Thus the action of this kernel on the initial, exponentiated data becomes 

/ au(x) 1 \ 2 / t=> u ( x / cos kt) a/ cos kt i \ 

ii(x,t,a) = K(t)(- = e * si "( 2fc *)oi^ [* i . (41) 



As in |10| . series expansion in powers of a and evaluation of l|41|l at a — produces a time series solution for u{x,t), 
since lim a ^ t, a) = hm a ^ - [ e au( - x ^~ akxta - n( - kt ' ) — l) = u(x,t) - fcxtan(fct). The series so obtained explicitly 
exhibits any difference between u and the particular solution incorporated into the exponential. 

1 00 1 / 1 \ ^ 

u(x,t) = kxt&nikt) H -V- -rf-smfct) -r~ (u (^r)) 1+J . (42) 

v ' y v 7 cos/ci ^ (1 + j)\ \k J dxJ v ^osfct;; v v 

Once having found these results so systematically - although perhaps somewhat laboriously - their interpretation 
is almost as simple as for the case of a constant driving term. As the reader may check by direct substitution, if 
u(x,t) is any solution of J3Jl, then ct J fct u ( co ^ fct , ^ tan kt) + kxt&nkt is a solution of © resulting from the same 
initial value data HDH2TJIH1. 
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The preceding modification of the propagator approach to can be employed for general pressure-gradients. If 
Up (x,t) is a particular solution to J2J), even if the pressure-gradient depends on both x and t, then 

11 (a;, t,a) = l f e »«(*.t)-««p(*.t) _ A (43) 

satisfies the linear equation 

d d 2 a \ , 4 a 



- (<9 x u p (x, *)) (\ + a^A ~ u p (x, t) U (x, t,a) = (44) 



if and only if du (x, t) jdt = u (x, t) du {x, t) / dx+g (x,t). Thus, evolution of any exponentiated initial data, expressed 
as exp (u (x) — u p [x]), can be effected by solving a linear equation for the kernel appropriate to and acting with 
that kernel on it {x, t = 0, a) . For certain problems, in addition to the g (x) — k 2 x case just discussed, this approach 
may be simpler to carry through than a direct attempt to construct the propagator for ifTjl. 

Breaking waves, shocks, and all that The formation of shocks under Euler-Monge evolution is well-known, 
and has a simple qualitative description in terms of characteristics (for example, see [141 lltjp . It is perhaps not so 
well-known that the series solution l|25[l leads to the same conclusions. In general, this may require some asymptotic 
approximations, but it is easy to find simple analytical examples. 

A particularly elementary illustration of a breaking kinematic wave is given by taking continuous, piecewise linear 
u profiles. First consider the case of the undriven Euler-Monge equation. The time power series sum ll'25ll for k = 
is easily evaluated in this case. Take one particular initial linear segment to be 

u (x) = a + Bx . (45) 

Then we have ^jrm -§^k u {x) n+1 = /3 n (a + fix), and so, for the time development of this segment, 

jr \ - a+l3x 

71=1 



u(x, t) = (a + fix) I 1 + J2 m n I = ■ (46) 



The segment becomes vertical at its "break-time" £break = 1//?, and while obvious, it is perhaps worth stressing that 
this is precisely the same as the radius of convergence of the time power series (|25|) . 

Continuously rising and falling combinations of such linear segments, forming a triangular "bump," with equal 
amounts of rise and fall, have conserved height and base length. Hence the triangle's area is constant. So too 
for the area under a continuous polygonal bump. This area conservation is just what you would expect from the 
Euler-Monge equation, written as a local conservation law: d t u(x,t) = d x (u 2 (x,t) /2). At least, this is what you 
would expect for continuous and differ •entiable functions. The fact that the same integrated conservation law also 
applies to continuous, piecewise differentiable functions may be surprising upon first encounter, but it is well-known 
to be true [CTITfi] . 

After the earliest such break-time, the continuous piecewise linear u profile will no longer be single- valued. If 
there are physical reasons to prohibit this, a standard procedure is to replace the continuous but multi-valued u 
with a shock discontinuity, positioned so as to preserve the previously mentioned area conservation (cf. the Rankine- 
Hugoniot rule, or the Lax entropy condition |2H|). 

When constant pressure-gradients drive the Euler-Monge equation, the result H25J) shows the same initial linear 
segment evolves as 

, a + 3 (x + \kt 2 ) , 

u (x,t) = kt+ — ~Y~pi — ' (47) 

with the same break-time as before, fbroak = !//?• When linear pressure-gradients are added to the Euler-Monge 
equation, IZfy shows the initial linear segment evolves as 

a + dx/coskt 

u{x, t) = kx tan kt H 3 , (48) 

cos kt — j sin kt 

with modified break-time given by ibrcak — \ arctan (k/ 8). Again, these break-times are just the radii of convergence 
for the respective time power series. 
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Another example is given by taking exponential initial data, for which we take the k — case to again see the 
wave evolution from the series Q25p. Suppose that initially u (x) = Ae x l L , then (x) n+1 — (^J^)" (Ae x ^ L ) n+ , 
so 

This series is familiar to anyone who has ever attempted to find by iteration the inverse function for the elementary 
equation x (y) — ye y , i.e. to obtain y (x) — xe~ y — x — x 2 + |x 3 — |a; 4 H — ■ • • = x~Y^_ x (-ra)" -1 jn\. Lambert 
and Euler studied the analytic properties of a generalization. There is a contemporary movement to standardize the 
function and name it after Lambert [24J . The resulting "Lambert W function" series definition, its log derivative, 
and its key functional relation are given by 

W{z) = -Y j ±n n - 1 {-z) n , z ±W{z) = ^P ( - W(z)exp(W(z)) = z. (50) 

n— 1 ' ^ ' 

This series converges absolutely for \z\ < 1/e. 

So the above solution Ij49(l for exponential initial data is 

u(x,t) = -^w(-^e*A . (51) 



t V L 

The radius of convergence of the series for u is therefore both t- and 2-dependent. Absolute convergence holds when 
^e x l L < 1/e . Actually, the Lambert function is well-behaved for all positive arguments, as is obvious from the 
inverse function, hence the exponential shape evolves backward in time without the formation of any singularities. 
But for forward time evolution, the wave always breaks, or forms a shock- front, and it does so at earlier times for 
larger x. 

The condition for the wave to break, for positive t, is just that the radius of convergence of the time power series 
is reached. That is fbreak = ^e~ 1 ~ x l L . At this time, for the given position, the slope of the wave profile becomes 
vertical, since 

r d , s u(x.t) 

L—u (x, t) = y \ - , (52) 

dx K ' l-tu(x,t) /L 1 v ' 

and W (—1/e) = —1 gives u (x,t = j^e~ 1 ~ x / L ) = L/t. The motion of the vertical face of the wave is simply given 

hy M eX /L = l 



or 



x(t) = L ^ln^ - 1 -hxt\ . (53) 

There is also the implicit definition of u (x, t) provided by when specialized to the case at hand: x + ut = G(u) = 
Lln(u/A), or 

u {x, exp(tu(x,t) /L) . (54) 

This functional relation is immediately seen to be the familiar one obeyed by the Lambert function, as in (|5U|) above. 
Indeed, our derivation of the series solution to the Euler-Monge equation for general data obeying (JSJ) may be thought 
of as a generalization of such series solutions to an infinite class of such functional relations. 

We plot below the A = 1 = L, k = case of u for various t. The exponential initial data u(x,0) = e x is 
shown as a black curve. With our sign conventions in the Euler-Monge equations, the wave moves to the left. For 
negative t the solution is also unbounded as x — > +oo, but non-singular for finite x, as shown by green curves in the 
right-most part of the Figure for various negative times. For positive t, vertical profiles occur where the curve has 
infinite slope, at the points indicated by small circles on the tips of the series-defined segments of u (x, t) shown by 
green curves in the left-most part of the Figure for various positive times. The series fails to converge for larger 
values of x, beyond these tips. (In more realistic situations, this may be indicative of some missing physics, such as 
would be represented by adding higher x derivatives in the original PDE, as in Burgers' equation, see |S].) However, 
the implicit definition of u (x, t) given by the functional relation JSJ serves to define upper branches of u (x, t) that 
are shown by the orange curves extending beyond the location of the vertical face. These upper branches may 
be thought of as the overhanging portions of the breaking wave. A pretty numerical example of a continuous, 
integrable, piecewise exponential wave is obtained by joining rising and falling exponentials, so that the initial data 
is exp (— |a;|). The time evolution of this example can be found graphically displayed on the website of one of the 
authors, (http://curtright.com/waves.html) 
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Evolution of the "Lambert front" for t = — 1/2 to +1/2 in At = 1/8 steps 

The location of the vertical face moves to the left as x (t) = —1 — lni, more slowly than the upper branch of the 
front, and reaches the origin at t = 1/e. (Recall the model is non-relativistic.) 

Moreover, we may extend the solution (|51|l of the Euler-Monge equation stemming from the exponential initial 
data, using the analysis of the preceding subsections, to find the constant body-force solution arising from the same 
initial data. 

u(x,t) = kt- - W ( - -^expi ( x + \kt 2 ] ] . (55) 



t \ L ' 1 L 
We also find the corresponding linear body-force solution. 

/ \ i i Lk ( At&nkt x \ , . 

u(x, t) — kx tan kt : — — W : exp — . (56) 

sin kt \ Lk L cos kt J 

Thus we have conditions for the wave profile to be vertical (signaling the local formation of shocks) whenever the 
arguments of the Lambert functions are equal to —1/e, i.e. 



''break 



L -r/r L -1 / 1, 2 \ , Lk ( -x \ 

= -7-e ' , tbroak = ~T ex P -T \ x + o fci broak , or tan «*break = ~T ex P 7 77 > ( 57 ) 

Ae Ae L \ 2 J Ae \L cos /cibreak / 

respectively, for the three situations considered. The latter two conditions are very implicit determinations of £break- 



5 Simple consequences of the extra dimension approach 

In this section a trivial method is given to construct an infinite number of solutions to the diffusion equation starting 
from a solution of the Monge equation, which might well be taken as one of the particular solutions of explicit form. 

This follows very simply from the existence of operators which commute with the diffusion operator. Consider 
the case where p = 0, i.e the free case. Then because Jj, ^ and the "boost" generator x-j^ — aJ^ all commute 

with the operator ^ — Q xda , we can find an infinite class of solutions of the diffusion equation based upon the 
solution u of the Monge equation by operating upon the "primitive" solution U(x, t, a) with arbitrary functions of 
these operators. (In fact only three such operators are required, as the diffusion operator is itself a combination of 
the others.) For example, consider 

G(^-)^-V(x, t, a) = G{u)^e a <^ , (58) 



9 



where G is an arbitrary function with no singularity at the origin, or 



G (|")| U (M,a) = G(u)^e™^ . (59) 
aa at at 



Note that in fact this latter class is really the same as the former, on account of the differential equation satisfied 
by u: All time derivatives may be replaced by derivatives with respect to x. In the case of generic g(x) ^ 0, the 
only operator which commutes with the diffusion operator is Jj. But there are exceptional cases. For example, 
when g(x) — k 2 x, the boost generator x-^ — a-^ also commutes. 

Suppose u is a particular solution of the Monge equation, in the case of the linear pressure-gradient term, where 
g(x) = k 2 x, such as one of the solutions given above: u = kxt&n{kt) or u — — kxcot(kt). Suppose further that 
another solution of the diffusion equation is sought of the form V(x, t, a) = v(x,t)e au{x ' t \ Then it may be verified 
that v(x,t) satisfies the equation 

If u is known explicitly then this equation can be integrated: In the case u = fcrtan(fci), v = — 1/ cos(fct), and in 
the case u = —kxcot(kt), v = 1/ sin(fct). 

6 Schrodinger and Euler-Monge relations 

The diffusion equation may be interpreted as a Schrodinger equation, upon complexification. To be explicit, consider 
the time-dependent Schrodinger equation for a harmonic potential in one dimension. 

di/>(x,t) d 2 ip{x,t) 1 2 2 

-^^ l ^x 2 —-V kx ^ t) - (61) 
We can rewrite this without change of content, together with its complex conjugate (with {x, a) both real) as 

dip(x + ia,t) d 2 , . . ,o , / n 

r = „ „ w(x + la, t) - -ik (x + la) wix + la, t) , 

ot axaa 4 

d%j){x - ia,t) d 2 - . .s,l.,2/ . ,2 1/ ■ .s / R0 ^ 

= - - ip[x — la, t) + —ik (x — la) ip{x — la, t) . (o2) 

at axaa 4 

Multiply the first by ${x — ia, t) and the second by ip(x + ia, t) and add to obtain 
d - d 2 - 

— {^p{x + ia, t)^p{x — ia, t)) — {}P{x + ia, t)ip(x — ia, £)) + k 2 ax (ip(x + ia, t)ip(x — ia, t)) . (63) 

Defining the complex point-split density p{x,t,a) — %jj{x + ia,t)ip(x — ia, t), we see that p obeys a diffusion type 
equation which is solvable in terms of a one-dimensional Euler-Monge equation with a linear pressure-gradient term, 
in the manner indicated above. Defining as previously x = \(x\ -I- ix%), a — \{x\ — 1x2), we find that p satisfies 
formally a two-dimensional diffusion equation with a complex harmonic potential: 

dt-dx\ + dx 2 4^i+z 2 )P- (64) 

This transformation of the Schrodinger equation is subtle, however, as now the arguments of the original wave 
functions are complex for real x\ and X2'- x Ha — ^ (1 ± i) [x\ ± 22). 

For real point-splitting, there are an infinite number of conserved charges and corresponding current densities 
for the free-particle Schrodinger equation, in an arbitrary number of spatial dimensions, as is well-known. These 
conservation laws follow immediately from 

d t ip (x, t) = iKV 2 ifj (x, t) , 9 f i/!(x,i) = -kV^(x,t) , (65) 

with k = For the real point-split probability density (just the density matrix in position-position representation), 

p (x, t, a) = 'ijj (x — a, f) ip (x + a, t) , (66) 

it follows that 

fitp(x,f,a)=V.J(x,t,a) ■ (67) 
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J (x, t, a) = intp (x - a, t) (V x ^ (x + a, t)) - in (V x ?/> (x - a, t)) i/j (x + a,t) . (68) 
However, as above, it is more in line with |l(Jj to write this current density as 

J (x, t, a) = in\7 a p (x, t, a) , (69) 

so that the infinite set of conservation laws become just the second-order equation 

d t p (x, t, a) = in (V x • V.) p (x, t, a) . (70) 

We recognize this as just the Fourier transform of the Wigner-Moyal equation |25l 1261 125] . With 

/(x,t,p) = J d n ae 2 ™ vp(x,t, a ) , (71) 

and H = up 2 , the free-particle Wigner-Moyal equation is given in terms of Groenewold's ^-product as 
-idtf (x, i, p) = -2ikv ■ V x / (x, t, p) = [H, f (x, t, p)] 4 =H*f-f*H, 

ff*/ = K(p-|iV x ) 2 /, /*^ = K (p+izV x ) 2 /- (72) 
Take one a gradient of (|7U[) to obtain 

d t J (x, i, a) = in (V x • V a ) J (x, t, a) . (73) 

Or take any number of a gradients. This gives a set of totally symmetric tensors conserved in precisely the same 
way. 

Tij...k (x,*,a) = d ai d aj ■ ■ •9 afc p(x,t,a) , d t Tij... k (x, i, a) = (V x • V a ) Tij... k (x, i, a) . (74) 

Of course, x derivatives also commute through d t — in (V x • V a ), so these may also be used to obtain conserved 
tensors, but as total x derivatives, they are not quite as interesting. On the other hand, the symmetric generators 

M ( y) = Xid Xj + xfi Xi - a.id a] - a^d^ (75) 

also commute with V x • V a , and may be used to construct other conserved densities. 

It is amusing that the conservation laws l|70|l are exactly of the form as in the introductory remarks of |10| , with 
complexification of the second-derivatives' coefficient. Unfortunately, if that coefficient is not purely imaginary, then 
the construction falters for real point-splittings. In other words, for real functions obeying the real diffusion equation 
this particular collection of conservation laws does not exist. Still, it is also amusing that one could indeed think 
of the point-splitting as introducing extra dimensions, and doing so by dimension-doubling. The wave functions 
involved in the currents are then just two distinct types of particular free wave functions in twice the number of spatial 
dimensions, namely -0(x, t, a) = (x + a, t) and ( x ? t, a) = (x — a, t). They each obey the extra-dimensional 
equations 

dtij) (x, t, a) = (V x + V a ) ^ (x, t, a) , 8$ (x, t, a) = -i^ (V x + V a ) ^ (x, t, a) . (76) 

There are other solutions of these free particle equations, of course, but as in ^U] we need only make use of particular 
solutions to encode the conservation laws. 



7 Generalized Bateman equation 

In previous investigations |29j . the Euler-Monge equation was used to obtain the prototypical Bateman equation 
|20]. This section addresses the question of what happens when the pressure term is added, and gives solutions of 
the resulting modified Bateman equation. 
First of all we notice that the equation 

d ( du\ (7?) 



dt 2 dx\ U dt 



holds whatever the pressure term is, provided it is dependent only upon x. This implies the existence of a scalar 
4>{x,t) such that 

du d(j> du d(j) 

W = fe' u m=Tf (78) 



ii 



That is, 

This and imply that </> satisfies a Bateman-like equation, in interaction with an external field: 

{(t>x) 2 <t>u - WxMtx + {(j>t?(t>xx = {(j>xfg{x) . (80) 

Note that the RHS preserves homogeneity in 0, but breaks homogeneity in t and x, in general. In the particular 
case where g(x) = k, a constant, and u — j {^~x + ^§-^, (f> is determined by l(75ll to be 4> — ^ (x + . Since 

any function of (f> will also satisfy the Bateman-like equation (as is clear from H79H ). we might as well take the square 
root of this for verification of (|80|l . 

If F((j>), G(tfi) are arbitrary differentiable functions of 0, then the general solution of the Bateman-like equation 
when the pressure-gradient is constant is given by solving implicitly for <j> the following equation: 

kt 2 

(x + — ) F{4>) + t G(4>) = constant . (81) 

Direct calculation establishes this result. We find 

—F , ktF + G , 



t - =. >Px , 



{x+^)F> + tG> ' F 



(82) 



and hence 



FG'-GF' 

p-2 — € > 

FG' - GF 1 



F 2 



K4>t - Hi . (83) 



Substitution of these last two expressions into the Bateman-like equation H80(l proves the claim. The implicit 
equation Ij81|l gives the general solution to (|8t)(l . for constant pressure-gradient, as it depends upon two independent 
arbitrary functions. 

A similar result can be found in the case of a linear pressure-gradient term, g(x) = k 2 x, which is the unique 
example that preserves homogeneity in x for H80|) . Here, the general solution is given by the implicit solution for 4> 
of the following equation: 

sinkt ,, coskt ,, . . 

F{6) -) GU) = constant . (84) 

x x 

Once again, direct calculation establishes this result. This is the second known case where a generalized form of the 
Bateman equation is integrable, the other being the so-called "two dimensional Born-Infeld equation" \M\ \\VZ\ . 

More general pressure-gradients appear to be less tractable for the Bateman-like equation, even implicitly, without 
resorting to the Euler-Monge form of the equation. Note that the recipe which worked to obtain implicit solutions 
of the driven Euler-Monge equation for arbitrary g (x) is not as useful for (|8Ufl . since thinking of t (x, <j>) instead of 
4> (x, i) does not linearize that equation. 



8 Summary 

We have considered here the effects of given pressure-gradients on the one-dimensional flows of the Euler-Monge 
equation, and on the corresponding driven diffusion/Schrodinger equations expressed in an extra dimension. The 
extra dimension method provides some additional insight and computational tools beyond those given by the method 
of characteristics and the technique of constructing implicit solutions, which we have discussed and compared. While 
we have touched on some elementary features of conservation laws for the free particle Schrodinger case, including 
their relation to the Wigner-Moyal equation, a full analysis of higher-dimensional driven flows along the lines of ^Oj 
remains to be given The effects of viscosity have also not been considered. (For a recent study including 

viscosity effects, see 8 .) We hope to return to these issues in a subsequent paper. 
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